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DifTusion induced by a thermal gradient in a granular binary mixture is analyzed in the context of 
the (inelastic) Enskog equation. Although the Enskog equation neglects velocity correlations among 
particles which are about to collide, it retains spatial correlations arising from volume exclusion 
effects and thus it is expected to apply to moderate densities. In the steady state with gradients 
c I ^ ' only along a given direction, a segregation criterion is obtained from the thermal diffusion factor A 

measuring the amount of segregation parallel to the thermal gradient. As expected, the sign of the 
5_j ■ factor A provides a criterion for the transition between the Brazil-nut effect (BNE) and the reverse 

Qh' Brazil-nut effect (RBNE) by varying the parameters of the mixture (masses, sizes, concentration, 

I solid volume fraction, and coefficients of restitution). The form of the phase diagrams for the 

BNE/RBNE transition is illustrated in detail for several systems, with special emphasis on the 
significant role played by the inelasticity of collisions. In particular, an effect already found in 
dilute gases (segregation in a binary mixture of identical masses and sizes but different coefficients 
' (— I ' , of restitution) is extended to dense systems. A comparison with recent computer simulation results 

^ ' shows a good qualitative agreement at the level of the thermal diffusion factor. The present analysis 

' generalizes to arbitrary concentration previous theoretical results derived in the tracer limit case. 
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; I. INTRODUCTION 

Controlling the mixing /demixing of granular media containing more than one species (polydisperse systems) is a 
problem faced by a wide range of industries. In some cases one would like to enhance the mixing effect while in other 
, situations it might be a desired and useful effect to separate particles of different types. Nevertheless, in spite of 
(~| its practical importance, the physical mechanisms involved in the segregation phenomenon are still not completely 
O ', understood [l| . This fact has motivated the development of accurate continuum models for polydisperse solid mixtures 
^ in order to offer a reliable description of the bulk behavior of these systems. 

' One of the most famous examples of (size) segregation in vertically vibrated mixtures is the Brazil-nut effect (BNE) , 
^f) where a relatively large particle (intruder) tends to climb to the top of the sample against gravity [H-Q . On the other 
J> ' hand, a series of experimental works [Gjj lli] have also observed the reverse buoyancy effect, namely, the intruder 
, can also sink to the bottom of the granular bed under certain conditions (the reverse Brazil-nut effect, RBNE). 
• Several mechanisms have been proposed to explain the transition BNE/RBNE, such as for example, void filling [2], 
[ convection [1, inertia @, and interstitial- fiuid effects @. Among the different competing mechanisms, thermal 
• • diffusion becomes the most relevant one when the granular system is vigorously shaken. Under those conditions, the 
motion of grains resembles the motion of atoms or molecules in an ordinary gas and so near-instantaneous binary 
\^ . collisions prevail. In this case, kinetic theory properly modified to account for the inelasticity of collisions may be 
' quite a useful tool to provide a reliable description of the kinetics and hydrodynamics of the system, and in particular 
• • ■ to analyze segregation in mixtures. 

Thermal diffusion (or thermophoresis in its single-particle manifestation [To| ) is the transport of matter caused by 
the presence of a thermal gradient. Due to the motion of the components of the mixture, a steady state can be reached 
; ^ ' in which the separation effect arising from thermal diffusion is balanced by the remixing effect of ordinary diffusion. 
I As a consequence, partial separation or segregation is observed and described by the so-called thermal diffusion factor. 
While this phenomenon has been widely studied in ordinary gases and liquids flTl, much less is known on thermal 
diffusion in the case of granular mixtures. It must be noted that for granular systems thermal diffusion can appear in 
vibrated systems even in the absence of an external imposed temperature gradient, as a consequence of inelasticity. 
In this case (energy supplied by vertical walls), the mean kinetic energy of the grains decays away from the source of 
energy giving rise to a (granular) temperature gradient. 

Previous theoretical attempts to describe thermal diffusion based on kinetic theory have been reported in the 
past few years. In the low-density regime, Serero et al. [l2l . [isj have studied the direct influence of inelasticity 
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alone on thermal diffusion segregation in the case of near-elastic particles |12| and finite degree of dissipation |13| . 
In particular, they find a novel effect, namely the fact that, even when the species differ only by their respective 
coefiicients of restitution a^, they may segregate when subject to a temperature gradient. However, they assume 
energy equipartition, which can o nly be considered as acceptable when ~ 1. In fact, the failure of energy 
equipartition in granular mixtures |14| has been widely confirmed by computer simulations [15] and observed in real 
experiments of agitated mixtures. Additional efforts for dilute granular mixtures have been made to assess the 
impact of the breakdown of energy equipartition on thermal diffusion [TtI . [18| . Interestingly, nonequipartition plays 
an important role since those results show that the relative position of the large particles 1 with respect to the small 
particles 2 is given by the sign of the control parameter (m2Ti/miT2) — 1, where and Tj are the mass and partial 
temperature of species i. While in an ordinary gas this sign is fixed only by the mass ratio (since Ti = T2), for a 
granular gas it also depends on the temperature ratio T1/T2 because of the lack of equipartition. This segregation 
criterion compares well with molecular dynamics simulations flTj . |l 9l| . 



In the case of dense granular mixtures, Arnarson and Willits [20| have determined the thermal diffusion factor for 
nearly elastic mixtures. However, their theory (which is based on the results of Jenkins and Mancini [2l|) differs from 
the theory for elastic particles f22'| only in the fact that it includes a sink term in the equation for the temperature 
and so no other inelastic effects are accounted for. Slightly different approaches [l^] based on kinetic theory have 
been invoked to get a segregation criterion j23 - [26| in the absence of a temperature gradient. In this latter case, the 
segregation dynamics of the intruder is only driven by the gravitational force. 

The purpose of this paper is to determine the thermal diffusion factor A of a moderately dense granular binary 
mixture described by the (inelastic) Enskog equation. Since the main interest here lies in the analysis of the effect of 
a thermal gradient on granular segregation, it will be assumed that no body forces (e.g., gravity) are present in the 
system. The segregation criterion is obtained from the factor A, which is explicitly given in terms of the parameters of 
the system (masses and sizes of particles, concentration, solid volume fraction, and coefficients of restitution). More 
specifically, the sign of A determines the tendency of the large particles to drift toward the cooler or warmer plate. It 
is apparent that the knowledge of the thermal diffusion factor allows one to analyze the origin of its sign and how it is 
related to the different parameters of the system. Previous theoretical results |27h29| | on thermal diffusion have been 
recently reported by the author of the present paper in the intruder limit case when the gas is driven by an external 
thermostat. The objective here is to extend the above results to arbitrary concentration and compare these theoretical 
results with some recent molecular dynamics simulations of a granular segregating binary system subjected to a 
temperature gradient. 

It must remarked that the factor A has been obtained from a solution [sil, |32| of the Enskog equation that goes 
beyond the quasielastic limit (and thus, it applies for a wide range of values of the coefficients of restitution) and takes 
into account the non-equipartition of kinetic energy. In this context, our theory subsumes all previous studies for both 
dilute [U [13 ^'^d dense [2^ [23l - [25j mixtures and, additionally assesses the influence of concentration on thermal 
diffusion without any restriction on the parameter space (for comprehensive review for mixture theories, see jssj). This 
is the main added value in this paper since our results can be relevant for comparison with experiments/simulations 
at finite densities. Moreover, it must be stated that the Navier-Stokes hydrodynamic equations are not actually 
solved in the present paper. Instead the results are conditional: if the temperature gradient has a given form, then 
the segregation criterion has a resultant form. For instance, in the case of a temperature inversion the results for 
segregation reverse. 

The plan of the paper is as follows. In Section [H] the thermal diffusion factor A is defined and evaluated by using a 
hydrodynamic description. The factor A is expressed in terms of the pressure^ and the transport coefficients , Un, 
and D12 associated with the mass fiux. All these coefficients have been explicitly determined from a Chapman-Enskog 
solution [31I, [S^l of the Enskog kinetic equation. The explicit forms of the pressure and the transport coefficients are 
displayed in Appendix [Xj The knowledge of p, L'f , Dn, and D12 allows one to determine the thermal diffusion A as 
a function of the parameter space of the system: the mass (TOi/m2) and diameter (cri/(T2) ratios, the concentration 
xi, the solid volume fraction 0, and the three independent coefficients of restitution of the binary mixture aij. 
In order to assess the impact of the different parameters on the segregation criterion, some special situations are 
analyzed and illustrated with detail in Section Hill In Section ITVl the form of the phase diagrams BNE/RBNE in 
the {ai/a2,mi/m2}-pl&iie is widely investigated by varying the parameters of the system in the case of a common 
coefficient of restitution (an = Q!22 — 0112 = a). Section IVl deals with the comparison between the Enskog theory 
and molecular dynamics results [30|] for the thermal diffusion factor. The paper is closed in Section IVII with a brief 
discussion of the results. 
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II. ENSKOG KINETIC THEORY FOR THERMAL DIFFUSION 

We consider a binary mixture of inelastic hard disks (d — 2) or spheres (d — 3) of masses nii and diameters 
(Tj {i = 1,2). Without loss of generality, we assume that ai > a2- The inelasticity of collisions among all pairs 
is characterized by three independent constant coefficients of restitution an, a22j and ai2 = Oi2i- For moderate 
densities, it is assumed that the velocity distribution functions /i(r,v;t) of each species are accurately described by 
the coupled set of inelastic Enskog equations [3^ [35j . Like the Boltzmann equation, the Enskog equation neglects 
velocity correlations among particles which are about to collide, but it takes into account the dominant spatial 
correlations due to excluded-volume effects. In the hydrodynamic description, it is assumed that the state of the 
mixture is characterized by the local number densities ni(r,t), the flow velocity U(r,t), and the (total) granular 
temperature T(r,t). These hydrodynamic fields are defined in terms of the velocity distribution functions fi as 

ni= J dvf,{v), pU = Y^mi J dvv/i(v), (1) 

"^ = Et / ^"^'^'(^)' (2) 

where p — mini is the total mass density, n — ni is the total number density, and V = v — U is the peculiar 
velocity. Assuming that there are no external forces acting on the mixture, the macroscopic balance equations for 
mass, momentum and energy can be derived from the Enskog equation. They are given by [3lj | 

Dtn, + n,V • U + = , (3) 

AU + p-iy-p-o, (4) 
Ar--y ^^ + -^(v-q + P: vu) = -cr. (5) 

n ^-^ nii an 

i 

Here, Dt = 9t + U • V is the material derivative, is the mass flux of species i, P is the pressure tensor, q is the heat 
flux, and ^ is the cooling rate associated with the energy dissipation in collisions. 

The constitutive equations for the irreversible fluxes ji, P, and q, and the cooling rate C have been recently obtained 
up to the Navier-Stokes (NS) order (first order in the spatial gradients) from the Chapman-Enskog solution [s^ to 



the Enskog equation [31|, |32[ . The results are 



.^^_!!!tlD,,Vlnn,-'^^l^D,2^lnn2-pD{VlnT, = -ji, (6) 



Pa/3 = pSap - V { VaC//3 + V^C/„ - • U (5„^ ) - • U 6al3 (7) 

q = ^T^Dq.iS/ In m - T^Dq^2^ In in - AVT, (8) 

C = 6"^ + C«v • U. (9) 

In these equations, Dij are the mutual diffusion coefficients, is the thermal diffusion coefficient, p is the pressure, r] 
is the shear viscosity, k, is the bulk viscosity, Dq,i are the Dufour coefficients, A is the thermal conductivity coefficient, 
is the zeroth-order cooling rate and Cu is a transport coefficient associated with first-order cooling rate. All the 
above quantities have been explicitly obtained by considering the leading terms in a Sonine polynomial expansion 
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A. Thermal diffusion factor 



As said in the Introduction, we are interested in analyzing segregation by thermal diffusion in a binary mixture. 
The amount of segregation parallel to the thermal gradient may be characterized by the thermal diffusion factor A. 
This quantity is defined in an inhomogeneous non-convecting (U = 0) steady state with zero mass flux (ji — 0) 
through the relation 

-A^ = ^.»f^V (10) 



dz dz \ n2 

where gradients only along the z axis (vertical direction) have been assumed for simplicity. Let us assume that the 
gas is enclosed between two plates where the bottom plate is hotter than the top plate, i.e., dz\nT < 0. In this 
geometry, according to Eq. (ITOl) , when A > the larger particles 1 tend to rise with respect to the smaller particles 2 
(i.e., dz \n{ni/n2) > 0). On the other hand, when A < 0, the larger particles fall with respect to the smaller particles 
(i.e., dz\n{ni/n2) < 0). The former situation will be referred here to as the Brazil-nut effect (BNE) while the latter 
will be called the reverse Brazil-nut effect (RBNE). 

We obtain now a relation of the type ^TU\\ from the balance equations. First, according to Eq. ©, the steady-state 
condition ji ^ — yields 

-{x^\xD\^^X2\2D\2)^Dl\ (11) 

where xi — nijn is the mole fraction or concentration of species i, 

dz In n, 

and we have introduced the reduced coefficients 

where i^q is an effective collision frequency defined below Eq. (|A5[) . Moreover, since U = 0, Eq. ([T]) clearly shows that 
the pressure tensor is diagonal for this state and so, Pa/3 = pSa^- In this case, the momentum balance equation Q 
reduces simply to 

The spatial dependence of the pressure p is through its dependence on the number densities and the temperature 
T. As a consequence, in reduced units, Eq. (1141) can be written more explicitly as 

- [xiPiXi + X2I32X2) = P* , (15) 

where p* — p/nT and 

Pi=T T—^P + ^+2^2-^, (16) 

ani x\ 0(p 0x1 



rr-idp ,^<f>2dp* dp* 

132 ^P + —-^ -xij—. (17) 

an2 X2 ocp 0x1 

Here, (f> = (f>i + 4>2 is the total solid volume fraction where 0i is the partial solid volume fraction of species i given by 

n,at (18) 



where T refers to Gamma function. 

The solution to the set of linear equations (fTTj) and (ITSt is 



p*Dl2~p2D^* p^Dj*-p*Dl^ 
' x,{(32Dl,-P,Dl2y ' X2{fi2Dl,-p,Dl2)' ^ ' 
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FIG. 2: (color online) Plot of the reduced coefficients (a) p(a)/p(l), (b) /3i(q)//Ji(1) and (c) /32{a)/ 132{1) as functions of the 
(common) coefficient of restitution a for hard spheres (d = 3) with xi = 1/2, axjai = 2, mi/m2 = 8 and a solid volume 
fraction = 0.1. 



Accordmg to Eq. (fTO|). A = A2 — Ai. Therefore, the thermal diffusion factor A is 

It is quite apparent that the influence of the parameters of the mixture on the sign of A is rather complicated, given 
the large number of parameters involved. For the sake of concreteness, we consider the region of the parameter space 
where P2D\i — /3i-Di2 7^ 0. In this case, the condition A = (which provides the criterion for the BNE/RBNE 
transition) implies [see the denominator of Eq. (|20p ] 

P*{xiDl^ + X2Dl^) = (xiA + X2p2)Dl*. (21) 

When the parameter set yields /?2-Dn - P1DI2 = 0, |A| 00 and BNE (RENE) appears if A > (A < 0). In any 
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case, the condition — ,Si-Di2 7^ covers practically all the parameter space of the problem. 

Equation (PT|) gives the curves delineating the regimes between BNE and RBNE. To get the dependence of the 
thermal diffusion factor on the parameters of the mixture, the explicit form of the transport coefficients and the 
equation of state is needed. These forms were evaluated in Refs. Isil . [s^ and some corrections to the expressions 
presented in these references have been done more recently in Ref . [37[ . The final correct expressions are displayed in 
Appendix \X\ for the sake of completeness. 

Figures [T] and [5] show the dependence of the diffusion transport coefficients {Du, D12, D]^), the pressure p, and its 
derivatives Pi with respect to the partial densities as functions of the (common) coefficient of restitution a = aij . We 
have considered here an equimolar mixture {xi = 1/2) with a solid volume fraction = 0.1 and mechanical parameters 
fi/c2 = 2 and mi/m2 — 8. To show more clearly the influence of inelasticity in collisions on mass transport and 
equation of state, all the quantities have been normalized with respect to their values in the elastic limit. We observe 
that the effect of collisional dissipation is in general significant, especially in the case of the diffusion coefficients. 



III. SOME SPECIAL LIMIT SITUATIONS 



The explicit form of the thermal diffusion factor A can be obtained when one substitutes Eqs. (jAl[) ~ (IA3l) for Df*, 
Z?*j^, and DI2, respectively, and Eq. (|A6p for p* (and its corresponding derivatives /?;) into Eq. (I^Ul) . This gives the 
dependence of A on the parameter space of the problem (mass and size ratios, mole fraction, solid volume fraction 
and coefficients of restitution). It is apparent that this dependence is in general quite complex. Thus, in order to show 
more clearly the different competing mechanisms appearing in the segregation phenomenon, it is first convenient to 
consider some special situations where a more simplified criterion can be obtained. 



A. Mechanically equivalent particles 

This is quite a trivial case since the system is in fact monodisperse {mi = m2, <ti — 02, an — 022 = 0112). In this 
limit case, Eq. (jAip shows that the thermal diffusion coefficient vanishes {Df* — 0), while D^i and DI2 are given by 
Eqs. (|A19p . Consequently, since the combination xiD^i + X2DI2 = 0, the factor A vanishes [see Eq. (PUI) ] and the 
condition (PT|) holds for any value of the coefficient of restitution and volume fraction. In this case, as expected, no 
segregation is possible. 



B. Dilute binary mixtures 

Let us consider a binary mixture in the low-density regime ((/)—>■ or, equivalently, riiaf — >■ 0). In this regime, 

p* — Pi — \, and 

ni^— + n2^— = C^"\ ni-^+n2^ =p. 22 
oni on2 oni on2 

Taking into account these identities, it is easy to get the explicit form of the transport coefficients from Eqs. (jAip - (jA3p . 
They can be written as 

Df* ^ i'^h - Cr' {^^111 - (23) 



1 \ / 



xiDl,+X2Dl,= \^iy*^^-Cj l^CDi* +xiji-^j , (24) 

where 71 = Ti/T (the partial temperatures Tj are in general different from the total temperature T) and v'^ and C* 
are defined by Eqs. (IA4p and (jASp . respectively, with Xij = 1- According to the above expressions, the criterion (PT|) 
becomes simply 



X1X2MC 



{v}, - C*)ix2 + xa){x2 + XiM) 
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FIG. 3: (color online) BNE/RBNE phase diagram for inelastic hard spheres {d — 3) with a = a\2 = ol22 ~ 0.8 in the tracer 
limit case {xi = 0) for two different values of the solid volume fraction 0. The solid lines correspond to the values derived from 
the relation (|3ip . while the dashed line refers to the results obtained when the gas is driven by an external thermostat, Eq. 
(|B12|) . Points above the curve correspond to A > (BNE) while points below the curve correspond to A < (RBNE). 



where M = mi/m2 is the mass ratio and 7 — T1/T2 is the temperature ratio. Since in general > C*? the solution 
to Eq. is simply 

1712 J2 

Although the explicit form of A derived here for a (unforced) dilute mixture differs from the one obtained when the 
mixture is driven (heated) by means of a stochastic external force the segregation criterion ([26l) (based on the 
sign of A) is the same as the one found in Ref. [31 . Note that if one assumes energy equipartition (Ti = T2), then 
segregation is only predicted for particles that differ in mass, no matter what their diameters may be. It must be 
remarked that the condition p6p is rather complicated since it involves all the parameter space of the system. As 
said in the Introduction, the criterion ([26]) compares well with molecular dynamics simulations (l7l | carried out in the 
tracer limit case {xi — )• 0) and is also able to explain [31 some of the molecular dynamics segregation results [31 
observed in agitated mixtures constituted by particles of the same mass density and equal volumes of large and small 
particles. 



C. Tracer limit in a dense binary mixture 

Let us consider now a dense binary mixture where one of the components is present in tracer concentration (xi — >■ 0). 
This problem is formally equivalent to studying the dynamics of an intruder immersed in a granular gas. The tracer 
limit case simplifies significantly the evaluation of the transport coefficients since, for instance, the dependence of 
the temperature ratio 7 = T1/T2 on the partial densities is only through the volume fraction [see Eq. (|B3|) ]. The 
explicit expressions for the diffusion coefficients in the tracer limit are given by Eqs. (IB4[) ~ (jB6l) . According to these 
expressions, Df* and D^i ^-''^ proportional to the concentration xi and so Eq. (|20| for the thermal diffusion factor A 
becomes 

^^^ f3x^'Dj*-p*{Dl,+x^'Dt2) ^ ^27) 



where /3 = p* 



p* = l + 2''-\22H^ + a22) (28) 
is the (reduced) pressure of the excess component and 



2<i-^dT{d/2) 



n2CT2 (29) 
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is the total solid volume fraction. Since (3 and D^^ are positive in the tracer limit, then the condition A = leads to 
the segregation criterion 



This criterion can be written more explicitly when one takes into account Eqs. (jB4l) - (jB6p with the result 

1 „ . sir , (i+w)'* M 



(30) 



1 + M' 



97 1 'y + M 6 / du 
j-M(3 + (j)-j- + -' ^ I ^ 



21 + MT 



It) + 



(31) 



where w = (Ti/(T2 is the size ratio and p,i is the chemical potential of the tracer particles [given by Eq. (IA12[) for c? = 2 
and Eq. (|A14p for d = 3]. It is apparent that in spite of the tracer limit case considered, the segregation criterion ((3T|) 
is quite intricate and so it is not easy to disentangle the impact of each effect (nonequipartition, dissipation, density, 
and/or mass and size ratios) on thermal diffusion. 

Thermal diffusion segregation of an intruder in a granular dense gas has been recently studied [13, HI] ■ In order to 
maintain the granular medium in a fluidized state, particles of the gas were assumed to be heated by a stochastic- 
driving force which mimics a thermal bath. This kind of forcing, which has been shown to be relevant for some two- 
dimensional experimental configurations with a rough vibrating piston (ssj , has been used in the past by many authors 
[ssij to analyze different problems, including segregation in granular mixtures [25| . Although previous experiments in 
vibrated granular mixtures [l6j have shown a less significant dependence of the temperature ratio Ti /T2 on inelasticity 
than the one obtained [1^ in systems heated by an external thermostat, some results (see for instance. Fig. 2 of Ref. 
[29']) derived for T1/T2 from this stochastic driving method compare well with molecular dynamics simulations of 
shaken mixtures [l9i] . This agreement suggests that this stochastic thermostat can be seen as a plausible approximation 
to modelize the experiments carried out in driven systems. On the other hand, more comparisons between results 
derived for driven and heated systems are needed before quantitative conclusions can be drawn on the reliability of the 
segregation conditions obtained from the transport coefficients derived with [29] and/or without [31, 32] an external 
thermostat. 

As expected, the external thermostat does not play a neutral role in the transport properties of the system [4Tj and, 
consequently, the criterion (|3ip differs from the one obtained in the driven heated case, Eq. (jB12l) . To illustrate more 
clearly these differences, a phase diagram delineating the regimes between BNE (A > when dzT < 0) and RBNE 
(A < when dzT < 0) in the {(Ti/(T2, m.i/TO2}-plane is shown in Fig. [3] for a = Q!22 — ot\2 = 0.8 and two different 
values of the solid volume fraction </). Significant quantitative discrepancies between the predictions obtained with and 
without a thermostat appear for small densities (0 = 0.1), since the effect of the thermostat being to reduce the size 
of the BNE region. Much less influence appears as the density of the gas increases since both results are practically 
indistinguishable at higher densities ((/) — 0.5). With respect to the influence of the solid volume fraction, we observe 
that in general the role played by the density of the gas is quite important since the range of size and mass ratios for 
which the RBNE exists increases with decreasing 0. 



D. Inelasticity-driven segregation 

In a previous theoretical approach [l3] based on a solution of the Boltzmann equation for nearly elastic particles, 
it has been found that segregation is induced by inelasticity. In other words, there is a separation between both 
species when they differ only by their respective coefficients of restitution. The authors explain the phenomenon as 
a consequence of the temperature gradient induced in the system by inelastic collisions, and relate the concentration 
gradient with the temperature gradient. These results have been subsequently extended to arbitrary degree of inelas- 
ticity [13'] . The above novel effect has been also confirmed more recently [12, by molecular dynamics simulations 
of a two-dimensional binary mixture kept fluidized by a vibrating base. 

In order to study the (pure) effect of inelasticity on thermal diffusion segregation, we consider the case m\ = mi 
and (Ti — (72 but different coefficients of restitution aij. Clearly, when all the coefficients of restitution are equal 
(ofii = a22 — OLYi), the system is monodisperse and so there is no segregation (A = 0). Figure S] presents plots of the 
thermal diffusion factor A as a function of the concentration x\ for different values of the coefficients of restitution at 
a density ^ = 0.2. As in the case of dilute gases (^ = 0) [l2], segregation in the presence of a temperature gradient 
can then occur due to inelasticity alone. In particular, segregation occurs even if one type of collisions is elastic (case 
(d)). While in this latter case the larger particles tend to accumulate in the warmer region, there is a change in the 
sign of A at a given value of the concentration in the other cases analyzed. In particular, the larger species tend 
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FIG. 4; (color online) Plot of the thermal difTusion factor A as a function of the concentration xi for a solid volume fraction 
(j) = 0.2 and different values of the coefficients of restitution: (a) an = Q22 ~ 0.5, 012 = 0.9; (b) an — 0.8, Q22 ~ 0.9, ai2 ~ 0.7; 
(c) an = 0.9, 022 = 0.8, ai2 = 0.7; and (d) an = 1, a22 = 0.5, ai2 = 0.75. Here, mi = m2 and tJi = a2. 
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FIG. 5: (color online) BNE/RBNE phase diagram for inelastic hard spheres [d = 3) with xi — 1/2, = 0.25, and two values 
of the (common) coefficient of restitution a = Uij . Points above the curve correspond to A > (BNE) while points below the 
curve correspond to A < (RENE). The solid lines are the results derived from Eq. (|31|l while the dashed line is the result 
obtained from Eq. (|3ip for a = 0.7 but assuming energy equipartition (Ti = T2). 

to move towards the colder plate when they experience more inelastic collisions than the other ones (an < a22)- 
Moreover, although not shown in the figure, our results also indicate a very weak influence of the volume fraction 
on the segregation process in this special case. 

IV. PHASE DIAGRAMS FOR THE BNE/RBNE TRANSITION 

Beyond the special limit situations considered in the previous Section, the thermal diffusion factor A [or, equivalently, 
the segregation criterion (j2ip ] depends in general on the following dimensionless parameters: the mass ratio mi/m2, 
the diameter ratio (Ti/(T2, the concentration xi, the overall volume fraction 0, and the coefficients of restitution an, 
a22 and ai2. For purposes of simplicity, henceforth the coefficients of restitution will be assumed to be the same 
for all combinations of collisions (i.e., an = 022 = ai2 = a). Moreover, we only consider the physical case of hard 
spheres (d = 3). This reduces the parameter space to five parameters. 

Next, we illustrate the form of the phase diagrams delineating the regimes between BNE and RBNE in the 
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FIG. 6: (color online) BNE/RBNE phase diagram for inelastic hard spheres (d = 3) with (/> = 0.25, ct = aij — 0.8, and three 



different values of the concentration xi: (a) xi — 0.1, (b) xi = 0.4, and (c) xi = 0.7 
A > (BNE) while points below the curve correspond to A < (RBNE) . 
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FIG. 7: (color online) BNE/RBNE phase diagram for an equimolar mixture {x-i — 0.5) of inelastic hard spheres (d = 3) with 
a = aij = 0.7 and three different values of the solid volume fraction 0: (a) — 0.1, (b) — 0.2, and (c) — 0.4 . Points above 
the curve correspond to A > (BNE) while points below the curve correspond to A < (RBNE). 



{CTi/tT2, TOi/m2}-plane as functions of the concentration xi, the overall volume fraction and the (common) co- 
efficient of restitution a. First, Fig. [5] shows a phase diagram for an equimolar mixture (xi = 0.5) at = 0.25 
(moderate density). Two different values of a have been considered {a = 0.9 and 0.7). In contrast to what happens in 
the intruder limit case (see, for instance, Fig. 4 of Ref. [1^), it is apparent that the main effect of coUisional dissipation 
is to reduce the size of the RBNE region. We observe that in general the RBNE (large particles tend to move towards 
the hot regions) is dominant for both small mass ratio and/or large size ratio. Moreover, in order to gauge the impact 
of the non-equipartition of granular energy on segregation, we have also included in Fig. [5] the phase diagram for 
a = 0.7 obtained from the segregation criterion ([21]) but assuming energy equipartition (Ti — T2). This has been a 
usual simplification in many previous theoretical works on thermal diffusion in nearly elastic systems (20l . [23 | . The 
comparison indicates a good qualitative agreement (at least in the region shown in the figure) between both results 
for not too large size ratios. On the other hand, quantitative discrepancies appear as the size ratio increases. In 
particular, although not shown in the phase diagram, when Ti = T2 the mass ratio becomes a two-valued function of 
the size ratio for values of (J\I<J2 ^ 2.6. This means that there exists a threshold value of a-ija^ above which no BNE 



11 



is predicted. It must be remarked that the significant influence of the different partial temperatures Ti on thermal 
diffusion found here is consistent with the molecular dynamics findings of Galvin et al. fsoj. These authors showed 
that non-equipartition driving forces for segregation are comparable to other driving forces for systems displaying 
comparable level of non-equipartition. Regarding this point it must be remarked that for systems where segregation 
is mainly driven by gravity (molecular dynamics simulations of Ref. (soj were carried out in the absence of gravity) , 
previous theoretical results |23 . [26l . [29j have clearly shown that non-equipartition has a weaker influence on segre- 
gation for thermalized systems (i.e., when d^T — )■ 0) than in the opposite limit (absence of gravity). This behavior 
qualitatively agrees with the experiments carried out by Schroter et al. 19]. 

Let us now analyze the effect of the concentration xi of the large particles on segregation. This is one of the main 
added values of the present paper with respect to previous studies focused on the tracer limit case {xi — >■ 0). Figure [S] 
shows a phase diagram for a = 0.8, </> = 0.25, and three different values of the mole fraction xi. We observe that the 
concentration of the mixture has significant effects in reducing the BNE region as xi increases. In particular, for a 
given value of the concentration, the transition from BNE to RBNE may occur following two paths: i) along constant 
mass ratio mi/m2 with increasing size ratio ai/<T2, and ii) along constant size ratio cti/(T2 with decreasing mass 
ratio mi/m2. Finally, Fig. [7] illustrates the influence of the volume fraction on the phase diagram for an equimolar 
mixture (xi = 0.5) at a moderate level of coUisional dissipation (a — 0.7). It is apparent that the role played by 
the density is quite important, especially for large size ratios. Surprisingly, and in contrast to the intruder limit case 
(see for instance Fig. 5 of Ref. 1231), the dependence on the solid fraction (j) is not monotonic: while the range of size 
and mass ratios for which the RBNE exists decreases with decreasing at moderate densities (lines corresponding 
to ^ = 0.1 and 0.2), the opposite happens at higher densities. Thus, at a given value of the size ratio, one has 
M"(0.1) < W'iQA) < M'''{0.2), where M^'icf)) denotes the critical mass ratio where the transition from RBNE to 
BNE occurs at density (j). 



V. COMPARISON WITH MOLECULAR DYNAMICS SIMULATIONS 



To the best of our knowledge, one of the few molecular dynamics simulations in which thermal diffusion has been 
isolated from the remaining segregation mechanisms has been reported by Galvin et al. [s^l . They consider a binary 
granular mixture constituted by frictionless inelastic spheres subject to an imposed temperature gradient. As expected, 
their results show in general segregation of particles according to their size and/or density (species segregation). Since 
these authors are mainly interested in assessing the role played by the non-equipartition of energy on segregation, no 
external forces like gravity are present in their simulations and hence the system is characterized by zero mean fiow. 

In the steady state, the granular temperature gradient between walls drives the segregation process. As in other 
experiments 0,0, the temperature profile demonstrates nonlinear behavior and exhibits a global minimum near the 
cold wall. In addition, Galvin et al. |30j also examine the profiles of solid volume fraction (/)(2) and partial densities 
ni{z) across the system to provide a quantitative measure of segregation (see, for instance. Fig. 7 of Ref. 

Although a direct comparison between the theoretical results derived here with those obtained in Ref. 



30 



)• 



301 for the 



partial density profiles ni{z) would require the numerical solution of the condition ji^z = along with the energy 
balance equation (O (in its steady state version), we will restrict our comparison to the thermal diffusion factor A. In 
fact, this quantity provides a more qualitative property of rii since its sign gives the tendency of each species to move 
upwards or downwards. In order to make a close comparison between theory and simulation for the thermal diffusion, 
let us consider the simulation data reported in Fig. 7 of Ref. (soj for a binary mixture with mechanical parameters 
o'i/cr2 = 2, mi/m2 = 16 and a (common) coefficient of restitution a = 0.9. According to the results displayed in 
this figure, it is quite apparent that A is a nonuniform function since it depends on z through its dependence on the 
volume fraction 4>{z) and the concentration xi{z) — ni{z)/n{z). Therefore, one can determine A((/)(z), xi(z)) from 
Eq. (pOt by using the local values of (f>{z) and xi{z) provided by the simulation as input parameters. This will give 
us the theoretical prediction of thermal diffusion across the system. The corresponding (local) value for A predicted 
by the simulations can be estimated from the temperature and species density profiles by numerically computing the 
derivatives dz\ii{ni/n2) and dzlnT at each point of the system. In this case, according to Eq. (jlOp . the value Amd 
given by the simulation is 

dzln (ni/n2) |md .„„^ 
Amd r. 1 ^1 , (32 J 

where the subscript MD means that these derivatives are obtained from the simulation data. 

The theoretical and molecular dynamics simulation results for A are shown in Table H] Regarding the sign of A, 
the comparison between theory and simulation shows a good agreement since both predict positive values for A in all 
the range of parameters {(f), xi) analyzed. Consequently, the more massive particles segregate preferentially toward 
the cool region, in qualitative agreement with the snapshot shown in Fig. 4 of Ref. [30| . At a more quantitative level. 
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4> 




Amd 


A 


0.33 


0.43 


1.09 


0.60 


0.45 


0.63 


1.21 


1.25 


0.50 


0.70 


1.02 


1.71 


0.51 


0.71 


0.83 


1.76 


0.48 


0.68 


0.90 


1.49 


0.41 


0.58 


1.03 


0.96 


0.30 


0.39 


1.26 


0.54 


0.20 


0.21 


1.80 


0.47 


0.10 


0.07 


2.44 


1.22 


0.08 


0.04 


2.29 


1.79 


0.06 


0.03 


2.19 


2.39 



TABLE I: Thermal diffusion factor as obtained from molecular dynamics simulations (Amd) and the Enskog theory (A) for 
different values of volume fraction (j) and concentration x\ for a granular binary mixture constituted by spheres with axja^ = 2, 
m\/m2 — 16, and a. = 0.9. Molecular dynamics results have been obtained by Galvin et al. [30l |. 

although theory and simulation compare well in the case of points near the minimum of temperature, there are in 
general discrepancies between theory and simulation. This quantitative disagreement can be due to the fact that while 
the expression (pn|) for A has been obtained up to the Navier-Stokes order (first order in the spatial gradients), the 
molecular dynamics simulations carried out in Ref. [soj clearly show strong gradients in solid fraction. Presumably, the 
numerical solution (beyond the Navier-Stokes domain) of the Enskog equation via the Direct Simulation Monte Carlo 
method (45| would give a better quantitative agreement with molecular dynamics simulations than the Navier-Stokes 
results reported here. This is a quite interesting problem to be addressed in the near future. 

VI. SUMMARY AND DISCUSSION 

The understanding of particle segregation within polydisperse, rapid granular fiows is still not well understood. 
The reason is perhaps twofold: first, there is a large number of relevant parameters involved in the description of 
the granular mixture; and second, there is a wide array of complexities that arise during the derivation of kinetic 
theory models. As previously mentioned, the two most common simplifications used in previous theoretical works 
on segregation have been to consider systems constituted by nearly elastic particles and an equipartition of granular 
energy. This paper has addressed the problem of segregation by thermal diffusion in a binary granular mixture 
described by the inelastic Enskog equation. The analysis has been based on a solution of the Enskog equation that 
covers some of the aspects not accounted for in previous studies. Specifically, i) it takes into account the nonlinear 
dependence of the transport coefficients on collisional dissipation (and thus the theory is expected to be applicable 
for a wide range of coefficients of restitution); ii) it considers the influence of the nonequipartition of granular energy 
on segregation; and iii) it applies to moderate densities. Consequently, the theory subsumes all previous analysis for 
both dilute [ll, [3 [13) 13 and dense ^2(1, .23- 25] gases, which are recovered in the appropriate limits. The results 
presented here generalize to arbitrary concentration previous results p7l - [29l | obtained in the tracer limit {xi — > 0). 

Among the different mechanisms involved in segregation, thermal diffusion (segregation induced by a temperature 
gradient) becomes the most relevant one when the granular system behaves like a granular gas. In the steady state 
with gradients only along a given direction, the sign of the thermal diffusion factor A [defined by Eq. ((TO)) ] provides 
information on the tendency of each species to move towards the colder or hotter plate. In this paper, the factor 
A has been evaluated by following two complementary approaches. First, by using the momentum balance equation 
(in the absence of gravity) along with the constitutive equation ^ for the mass flux, A has been expressed [see 
Eq. pO|] in terms of the pressure p (and its derivatives with respect to the partial densities Ui) and the transport 
coefficients Dn , D12, and Df associated with the mass flux. Then, the forms of the pressure and the diffusion 
transport coefficients have been explicitly obtained from a Chapman-Enskog solution of the Enskog equation [sTI, . 
This finally gives A as a function of the mass and size ratios, the concentration, the solid volume fraction, and the 
coefficients of restitution. In particular, the condition A = [see Eq. ((2T|) ] provides the segregation criterion for the 
transition BNE<^ RBNE. 

In general, the segregation criterion (|2ip presents a complex dependence on the parameter space of the system. 
In order to disentangle the impact of the different parameters on thermal diffusion segregation, some special cases 
(dilute gas, tracer limit, . . .) have been separately studied. An interesting new case corresponds to the segregation of 
a binary mixture of granul ar p articles that differ only by their coefficients of restitution. This novel effect was first 
predicted by Serero et al. |l2l . [l3| from the Boltzmann equation (low-density gas) and has been recently confirmed 
by molecular dynamics simulations of hard disks {43] . The results obtained here for dense granular binary mixtures 
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confirm also the existence of segregation induced by an inelasticity difi'erence (see Fig. |4]). Moreover, our results also 
show in general a weak influence of the volume fraction on thermal diffusion for this special situation. 

A systematic study of the form of the phase diagrams BNE/RBNE in the mass and size ratio plane has been carried 
out in Section HV] for hard spheres in the case = a. Regarding the influence of collisional dissipation on the form of 
the phase diagrams, the results indicate that the influence of a on A is quite important, the main effect of dissipation 
being to increase the size of the BNE region (see Fig. [5]). In addition, we also conclude that the role played by the 
nonequipartition of granular energy on segregation is quite relevant, especially as the disparity of masses and/or sizes 
increases. This result is consistent with recent molecular dynamics simulations 30]. With respect to the influence 
of the concentration xi, our results show that in general the main effect of xi is to reduce the BNE region as the 
concentration of the large particles increases (see Fig.[S]). Finally, we also observe that the form of the phase diagrams 
changes significantly with the volume fraction specially at large size ratios (see Fig. [7]). 



By extending the intruder limit analysis [27|,|29| to arbitrary values of concentration xi, comparisons with molecular 
dynamics simulations become practical and this allows one to assess the reliability of the Enskog kinetic theory to 
characterize thermal diffusion segregation. To make some contact with the molecular dynamics results of Galvin et 
al. 'sd] , we have compared the kinetic theory predictions for the thermal diffusion factor A for different values of the 
concentration xi and volume fraction (f> (for the system a\lu2 = 2, mxjmi = 16 and a = 0.9) with those obtained 
from the simulation data by numerically evaluating the derivatives d z\"i^(n\ I n-i) and Q^lnT at the points of the 
system corresponding to the same values of x\ and (/). The comparison between theory and simulation shows a good 
qualitative agreement since both predict the same sign of A for the different points analyzed. In addition, at a more 
quantitative level, although molecular dynamics simulations show strong gradients in the bulk region (and so they 
go beyond the linear domain of the Navier-Stokes description) , the theory compares reasonably well with simulation, 
especially in the region close to the minimum of granular temperature. It is important to remark again that the 
quantitative discrepancies between theory and molecular dynamics simulations are essentially due to the limitations 
of the Navier-Stokes results rather than the assumptions inherent to the Enskog kinetic equation (molecular chaos 
hypothesis). 

Certainly, the derivation of kinetic theory models for segregation flows in polydisperse systems is perhaps one of the 
most important open challenges of granular gas research. The theoretical results reported in this paper cover part of 
this challenge, at least in the case of the thermal diffusion segregation. On the other hand, the present theory has some 
important restrictions. First, given that the Enskog equation still assumes uncorrelated particle velocities (molecular 
chaos hypothesis) , it is expected that the kinetic theory for thermal diffusion only applies to moderate densities (solid 
volume fraction t ypic ally smaller than or equal to 0.25). However, despite this limitation, there is substantial evidence 
in the literature |46j on the reliability of the Enskog theory to accurately describe macroscopic properties (such as 
transport coefficients) for a wide range of densities and/or collisional dissipation. Another important limitation is 
that the segregation criterion derived here has been obtained by using the first Sonine approximation for the diffusion 
transport coefficients. Recent results j2^] for the tracer limit clearly show that the accuracy of the first Sonine solution 
can be questionable for small values of the coefficients of restitution and/or disparate values of the mass and size 
ratios. The infiuence of the second Sonine correction to the transport coefficients is an interesting open problem to 
be carried out in the near future. This will allow us to offer a segregation theory that can be reliable even for extreme 
values of dissipation or mass and size ratios. 
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Appendix A: Expressions of the mass flux transport coefficients and pressure 



In this Appendix we provide the expressions of the (reduced) transport coefficients £'12, and Df* associated 
with the mass flux and the hydrostatic pressure p*. These quantities are involved in the evaluation of the thermal 
diffusion factor A. 
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The expressions of the reduced coefficients Df* , D^^ , and can be written as [31 
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In these equations, 7^ — Ti/T, CJ* 
mi /{mi + mj), and 



C'"°V^o, P* — p/{nT), Xij is the pair distribution function at contact, Mij 
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where 9i = miT/moTi and mg = (TOi+m2)/2. The partial temperatures Ti and T2 are determined from the condition 
Ci°^ — C2''' ~ C'''^\ where the expression of c}f^ is 
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where I'o = 't.o'ij ^ \J2T jma. The reduced pressure p* is given by [3l[ 
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The expHcit form of the transport coefficients D*i and requires the knowledge of the the quantities Inj . These 
parameters are given in terms of the functional derivative of the (local) pair distribution function Xij with respect 
to the (local) partial densities ni [see Eq. (Cll) of Ref. [3l|]. Given the mathematical difficulties involved in the 
determination of the above functional derivatives, for the sake of simplicitjjthe parameters Inj are chosen here to 
recover the results derived by Lopez de Haro et al. for elastic mixtures [l^ (see Appendix C of Ref. [s^])- The 
quantities Inj are the origin of the primary difference between the standard Enskog theory and the revised version 
for elastic collisions ^8|. They are zero ii i = i, but otherwise are not zero. These quantities are defined through the 
relation J32|] 
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where fii is the chemical potential of species i and B2 = /dT{d/2) [B2 = § for disks (d = 2) and B2 = ^ for 
spheres {d — 3)]. Since granular fluids lack a thermodynamic description, the concept of chemical potential could be 
questionable. As said before, the presence of /i^ in our theory is essentially due to the choice of the quantities luj. 
Given that the explicit form of the chemical potential must be known to evaluate the diffusion transport coefficients, 
for practical purposes, the expression considered here for /Ji is the same as the one obtained for an ordinary mixture 
of gases {aij =1). Although this evaluation requires the use of thermodynamic relations that only apply for elastic 
systems, we expect that this approximation could be reliable for not too strong values of dissipation. More comparisons 
with computer simulations are needed to support the above expectation. 

Taking into account Eq. (IA7I) . the nonzero parameters /121 and /122 appearing in Eqs. (jA2p and (jA3p are given by 
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Note that for mechanically equivalent particles (mi — m2, ci — 0-2), -^121 — I122 — 0, as expected since the standard 
and revised versions of the Enskog equation lead to the same Navier-Stokes transport coefficients for a monocomponent 
gas 

In the case of hard disks {d — 2), a good approximation for the pair distribution function Xij is [50| 
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where — niiraf /4 is the solid volume fraction for disks and 
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The expression of the chemical potential /i^ of the species i consistent with the approximation (lAlOp is [5l| 
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where Xi{T) is the (constant) de Broglie's thermal wavelength 1521 . In the case of hard spheres {d — 3), we take for 
the pair distribution function Xij the following approximation [53| 
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where — niTTaf/6 is the solid volume fraction for spheres. The chemical potential consistent with (jA13[) is [52 
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According to Eqs. (IAl[) - (jA3[) . the diffusion transport coefficients are given in terms of the derivatives of 7^ with 
respect to the partial densities Ui. In terms of the temperature ratio 7 = T1/T2, the partial temperatures 7^ are 
defined as 
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The dependence of the temperature ratio 7 on ni and n2 is through its dependence on the concentration xi and the 
volume fraction (j). As a consequence, 

where 0^ is defined by Eq. (jl8p . The derivatives 9^7 and 93;j7 can be obtained by taking the derivatives with respect 
to 4> and xi in the condition (^l = C2 ■ This yields the relations 
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This allows us to express the derivatives 9^7 and dx-^j in terms of the temperature ratio 7 and the parameters of the 
mixture. 

It is apparent that the (reduced) transport coefficients have a complex dependence on the mass and size ratios, the 
concentration, the volume fraction, and the coefficients of restitution. A simple but nontrivial case corresponds to a 
binary system constituted by mechanically equivalent particles (mi — 1112, (Ji =01, cm = 0122 = 0^12) ■ In this case, 
Xij = Xj 7j = 7 = I7 P* = 1 + 2'^~^x^(l + a), and the parameters =0. As a consequence, Eqs. (|Aip - (jA3l) simply 
reduce to D\* = and 

i?t2 = -^^n, 0\1={^*D-\C^ [l-a;i(p*+(^a^p*) + 2a;i(p*-l)(l + 09^1nx)], (A19) 

where 



, x(l + «), C = ^-^7T7i^x(l - (A20) 



Appendix B: Tracer limit case 

The explicit forms of the diffusion transport coefficients in the tracer limit case {xi — ^ 0) are displayed in this 
Appendix. In this limit, 72 = 1, 71 = 7(0) and p* = 1 + 2''~^X220(1 + 022). Moreover, since the dependence of C'"', 
p and 7 on the partial densities is only through the volume fraction 0, one gets the simple relations dmC^^"^ = dmP = 
9„i7 = 0, 



n2^^C> l + cj^^^ I (Bl) 



The explicit expressions for the transport coefficients in the tracer limit can be easily obtained from Eqs. (|Aip - (jA3P 
when one takes into account the identities (IBll)-(IB3l). The result is 
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Here, (3 = p* 



, /ii is the chemical potential of the tracer particles and 
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As said in Sec. Illli in order to maintain the granular medium in a fluidized state, previous works [27l[29| considered 
the presence of a stochastic external thermostat. The corresponding expressions for the transport coefficients are [H^l 
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Taking into account Eqs. (|B9[) - (jBll[) . the segregation criterion ([3D]) becomes 
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